Fluctuations in jet momentum as an energy loss probe 
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We investigate the measurement of the distribution of away-side jet pr as a way to probe the 
energy loss mechanism in heavy ion collisions, and constrain the properties of the medium created 
in heavy ion collisions. We define an observable related to the fluctuation of the energy deposited in 
the medium, and show that competing models of parton energy loss give differing values and scaling 
of this observable. We also argue that the scaling of this observable with system size can be related 
to the medium's partonic density. We then give a qualitative discussion of how the measurement 
we suggest can be performed and used to determine the parton energy loss mechanism and possibly 
the system's partonic density. 

PACS numbers: 13.87.-a, 12. 38. Aw, 25.75.-q, 24.60.-k 



Jet suppression has long been regarded as a promis- 
ing signature for deconfinement/ Quark Gluon Plasma 
(QGP) production. It is very reasonable to suppose that 
partonic energy loss in a deconfined medium is substan- 
tially greater than a medium where quarks and gluons 
are locked in color- neutral bags . 

Perturbative Quantum Chromodynamics (pQCD) can 
be used to analyze jet energy loss quantitatively, and 
to explore the features that suppression in a deconfined 
medium should exhibit HHHIil. 

The convincing observation of jet suppression in RHIC 
A-A collisions, and the determination, via d-Au control 
measurements, that in-medium properties rather than 
initial state effects are behind the suppression is a cor- 
nerstone for the evidence that RHIC produced a "new 
state of matter" H 0, S H . 

However, this consensus is undermined by the ambi- 
guity over the correct approach to describing parton en- 
ergy loss in terms of QCD. Several models, treating jet- 
medium interactions through mutually inconsistent ap- 
proximations, have been shown to successfully describe 
available experimental data 0, ^| . 

As a result, translating the jet energy loss observation 
into a quantitative determination of the properties of the 
medium created at RHIC, such as an estimate of the 
initial partonic density or temperature, has been prob- 
lematic. 

Part of the reason for this ambiguity is that traditional 
jet identification and analysis tools are inapplicable to the 
heavy ion enviroenment, since the background of hun- 
dreds of low energy tracks makes it impossible to distin- 
guish every "jet" track from background low-momentum 
tracks independent of the jet. 

This makes quantitative studies of jet broadening and 
jet structure in A-A collisions, capable of testing energy 
loss models, much more difficult, from a phcnomcnolog- 
ical and experimental point of view, than jet studies in 
e-e or p-p collisions. 

QCD-based jet energy loss descriptions can be roughly 
divided in two classes: 

i) Thin plasmas, where it is assumed that the system 
size is of the order of the parton's mean free path, and 



hence opacity (the mean number of hard collisions n) 
becomes a suitable expansion parameter. The dynam- 
ics within this approach can be extrapolated from an S- 
matrix like formalism. The energy loss can therefore be 
calculated from a finite number of Feynman diagrams de- 
scribing hard partonic collisions. The most highly cited 
application of this a ppr oach is known as GLV (named 
after the authors) llllllllil. 

ii) Continuous absorption, where it is assumed that the 
mean free path is much smaller than the path the jet has 
to traverse before hadronizing. Hence, energy loss rate 
can be represented as occurring over "many" collisions 
occurring over each element of the parton's path within 
the medium. The total energy loss can be calculated by 
integrating the resulting stochastic rate equation 0, 0|. 
This approach was independently developed in 0| (also 
named after the authors, BDMPS) and [3- Recently, it 
was extended to make direct contact with finite temper- 
ature QCD via the so-called "AMY" formalism 

Note that the opacity parameter within the GLV for- 
malism is related to the partonic density in BDMPS and 
the medium temperature of AMY, since the mean free 
path parameter Xg can be related to the length of the 
parton's in-medium path L, the interaction cross-section 
a and the partonic density n via 
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The density in an ideal gas of light particles is (for np 
light flavors and three colors) is then given by 



(2) 



Similarly, the color screening mass, appearing as a cut- 
off in GLV minimum gluon virtuality [l^l, is related to 
the medium temperature and the strong fine structure 

constant a. via 



2 = 4™,T2 
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The crucial difference between the two approaches, used 
to devise our experimental observable, is that in the thin 
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plasma approximation the probability of no interactions 
remains finite, even if higher order terms in opacity are 
considered. The continuous absorption approach, on the 
other hand, is based on a stochastic equation where the 
transition rates include infinite re-summations of dia- 
grams. The probability of no parton-medium interaction 
within this ansatz is manifestly zero by construction. 

It is apparent that the probability density function of 
the hard parton's momentum as a function of distance 
traveled behaves in a very different way in the two ap- 
proaches. 

Within the GLV ansatz, it maintains a primary peak, 
centered around the same pT and having the same width 
as immediately after the initial scattering. The ampli- 
tude of this primary peak decreases with time, and even- 
tually goes to zero when the length traversed is much 
larger than the mean free path (a regime where the "thin 
plasma" approximation breaks down), but it's position 
and width should not change beyond effects due to frag- 
mentation and primary jets reinteraction (common to all 
systems). As the interaction probability grows, a sec- 
ondary peak develops. Since the collisions within the 
GLV formalism are typically hard, this peak should be 
well removed, in phase space, from the primary peak. 

This picture is quantitatively investigated in Fig. ^ 
In the left panel, we plot the probability density func- 
tion for the surviving fraction of the initial momentum, 
calculated following the prescriptions in 
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where pini is the initial parton momentum and ^ is 
given, to first order in opacity, by 
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/O 7r(qi+M2)2 16x2pL + (k-q)ii2 
Here the A:^^^. is given by kinematic constraints 

and N is a normalization constant to ensure the proba- 
bility density function is correctly normalized 
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As can be seen, re-scattering manifests itself by a sec- 
ondary peak, whose area saturates at unity (probability 
of collision= 1) as the size of the system becomes sig- 
nificantly larger than the mean free path. Multiple in- 
teractions (beyond first order in opacity) will change the 
secondary peak (Eq. El but, due to the "thin plasma" 
approximation, can not affect the center and the width 
of the primary peak (the (5- function in Eq. . 



The right panel of Fig shows the same distribution 
where the integral in Eq. [3 is restricted to a small angle 
between the medium parton k and the virtual gluon q 
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and only the full angular distribution is normalized to 
unity. It shows that the secondary peak is mostly at 
a large angle w.r.t. the original jet, significantly larger 
than a sensible experimentally defined jet cone. 

Hence, if only hard particles within a jet-cone angle of 
^ 0.2 w.r.t. the original parton direction are measured, 
it should be the case that either no jet particles are found 
(the jet has been completely absorbed by the medium) 
or a jet that has approximately the original parton's mo- 
mentum. 

Within heavy ion collisions, of course, fragmentation, 
together with background tracks, makes the definition of 
"jet particle" experimentally ambiguous. Furthermore, 
if we only observe one jet, we do not know the initial 
momentum of the original parton. 

If, however, we observe a pair of azimuthally corre- 
lated jets, the original parton momentum can be inferred 
by the higher momentum "near-side" jet (probably close 
to the surface, due to the high opacity of the medium). 
Provided the parton hadronizes outside the medium, the 
broadening and hadronic momentum distribution due to 
fragmentation should be similar to that in d-Au colli- 
sions. 

Hence, the GLV approach predicts that, while the 
probability of having two azimuthally correlated jets de- 
creases, once they are found they also maintain the mo- 
mentum correlation expected from d-Au collisions. 

The picture in the BDMPS/AMY formalism if very 
different. The stochastic evolution ensures that, instead 
of a secondary peak developing, the primary peak broad- 
ens and eventually disappears as the jet parton loses it's 
high momentum and becomes part of the medium. The 
fiuctuation in the jet's momentum, therefore, is expected 
to increase with the distance traversed in the medium by 
the parton. 

Fig. 121 shows the fraction of the original momentum 
plotted against traversed length for a variety of medium 
temperatures and system sizes. The calculation was done 
using a rate equation |lCt IT^ . where the approximate 
solution is given by 

Pamy (p, t) K. 

deD'^ie,p,t) P%,y ip + e,0) + DS{e,p,t) {p + e,0){8) 

where the Kernel is related to the transition rates 

r(p, uj, t) calculated from finite-temperature QCD within 

the AMY formalism (ij 

r roc 1 pt 

L>«/^'(e,p,t) = exp -/ duj X dt' r{p,uj,t') 

I J-oc J Jo 

n/ ^'^^ / dt'T{p,uJ,t')5[e-Y,^^ 
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FIG. 1: Left panel: probability density function for the partonic p/pinuiai calculated in the GLV formalism, to first order in 
opacity, as a function of distance traveled by the parton . Right panel: The secondary (rescattered) peak has been restricted 
to be within an angle < dmax of the original trajectory. See text around Eq. Qfor more discussion. 



the initial condition is given by a well-defined initial 
quark momentum 

PAMviP^t^O) - S{p-p,m) (10) 

In addition, to prevent parton-medium ambiguity, we re- 
quired that the parton has an energy above a minimum 
initial px of 2 GeV. 

It can be seen that the initially sharp probability peak 
broadens and reaches a maximum width. 

Subsequently, the probability density function concen- 
trates around the minimum pT'. This peak is a conse- 
quence of our introduction of a minimum pT in the prob- 
ability density function. It is not expected to be physical 
but is rather an artifact of our "jet" definition. It's ap- 
pearance coincides with the parton becoming part of the 
medium, and the jet becoming unobservable. 

As gluon splittin g a. t small angles dominates in the 
AMY approach [l^, [ij, , the broadening seen in Fig. 
121 is independent of the jet-cone angle chosen, provided 
this angle is large enough to capture the whole jet and 
small enough to keep uncorrelated jets out. therefore, 
a jet-cone angle similar to that defined within the GLV 
approach will also work here. 

Thus, we have shown a qualitative difference between 
GLV and AMY energy loss Ansatzes. Within GLV, the 
jet is absorbed through a decrease of it's amplitude, as 
expected from an approach derived via an "S-Matrix" like 
framework. Within BDMPS/AMY, the jet disappears 
after an initial broadening, again as expected from the 
stochastic nature of the underlying ansatz. 

These considerations lead us to propose the following 
observable: We trigger on 

• An near-side hard particle with momentum Pxriear 

• A away-side hard particle with a smaller but still 

jet-like momentum Ptngger < Pxaway < Prnear, 



where 

TT - 69/2 < 0near " Oaway < TT + 69/2 

we compute total hard momentum in a cone with angle 
69, defined as 

P^>Ptr>gger,|ej-e„ear|<<5e 

Pror{pTrmn,69)= ^ Pt, (H) 

3 

we then calculate the equivalent away-side. 

pi.>Ptrigger,\9j-9au,ay\<Se 

PZT'iPTrmn,59)= ^ Pt, (12) 

3 

(note that the number of such particles j in each case 
may well be 1) 

Our observables are then the average and event-by- 
event fluctuation of z defined as 

z = p^rvPtoT (13) 

and it's dependence on system size (centrality and nu- 
cleus type). Since we are using Pp^T^ as a gauge for the 
initial momentum of the parton, we require that z < 1. 
Hence, if in an event z > 1 (an unlikely but possible out- 
come) this variable should be flipped {z — > z"^) before 
being included in the statistics. 

The considerations of the previous paragraphs show 
that the "thin plasma" approximation demands that 

(2)j-.4 ~ <-U-An (") 
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FIG. 2: probability density function of p/piniuai as a function of distance traveled by the parton with the AMY formalism, for 
a variety of temperatures (particle densities). All distributions are normalized to one, and require that pr > PTmin to avoid 
in-medium contamination. The unphysical peak that develops at late times at pTmin is a consequence of this requirement. 



(where {{AX)^) = (X^) - {xf) 

Furthermore, GLV predicts no dependence with coUi- 
sion system size of either the average or the fluctuation 
in z. A higher parton density/opacity will make collect- 
ing the statistical sample required for measuring these 
observables harder, but once such a sample is obtained 
the results should follow the scaling of Eqs. ll4land[T31at 
any finite opacity. 

If, on the other hand, 

^')A-A < Wj-.4. (m) 

((^=>').-. > (<^^''>.-.. <") 

and ^(Az)^^ exhibits a strong dependence on system 

size (with a maximum in system size similar to the one 
present in Fig. \^ than an AMY/BDMPS hke ansatz is 
more appropriate for describing jet energy loss within the 
system. In this case, the system size where the maximum 

(^(Az)^^ is observed can be related to the tempera- 

ture/partonic density of the created medium. 

The leading source of systematic error for this mea- 
surement is the presence of uncorrelated jets within the 
system. Hence, for this method to work, the important 
criterion in the definition of Ptrigger is to make sure that 
the probability of having uncorrelated hard particles in 
the same events with P > Ptrigger is small, so 

((Az)^) ^ ((zf ) (18) 

To investigate the suitable ptrigger and 56 for this mea- 
surement, we have analyzed the observables of interest 
using HIJING [13 - 



The result is shown in Fig. 13 for a variety of jet cone 
angles and trigger momenta. As can be seen, if the jet 
cone angle is < 0.6 radians and pTmin < 4.5 GeV, the 
and ^(Az)^^ converge independently of the 
precise jet cone angle value. 

The appearance of a strong jet cone angle dependence 
of (z),((Az)^) at PTmin = 4.5 GcV helps define an 

appropriate upper limit for pTmin (as per Eqs. 1111 and 
I12|l . At that PTmin , the probability density for z is dom- 
inated by the contamination of the chosen jet by inde- 
pendent hard processes (neighboring jets). Of course, HI- 
JING does not include the significant quenching observed 
at RHIC, so the pTmin limit can actually be pushed for- 
ward when analyzing experimental data. 

The onset of BDMPS/AMY regime should manifest it- 
self by the broadening of the z distribution. Such broad- 
ening should increase ( (Az)^ ) for at least a sub-set 

of system size bins, up to a maximum value, fixed by the 
definition of z to ^ 0.5. 

As Fig. [21shows, the system size at which {^(Az)^^ 

saturates to ^ 0.5 can be directly related to partonic den- 
sity (or temperature in the AMY approach). For T=540 
MeV it's ~ 2 fm, for T=360 MeV it is ~ 4 fm and for 
T = 180 MeV it is ~ 6 fm. Hence, measuring this sys- 
tem size experimentally would remove the density/size 
ambiguity inherent in the interpretation of average jet 
suppression measurements. 

Such a measurement should be experimentally feasible: 
At the desired system size (be it centrality or colliding 
nucleus type), angular correlation of the produced jets 
should be nearly as high as in p-p and d-Au systems, 
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FIG. 3: Left panel: (z) and ((Az)^) calculated with HIJING and plotted for a variety of jet cone angles 56 and trigger momenta 
PTmin (See eq. ll3lllll and ll2l for definition). Empty symbols refer to d-Au, while full symbols are Au-Au. Errors were calculated 
using the Jackknife method 18|. Right panel: A schematic representation of the expectation for these observables in the GLV 
and AMY Ansatzes. The dashed-line box in the bottom panel shows the system size of ((Az)^) saturation, related to the 
medium's partonic density 



Since the distance traversed by the parton is not enough 
for a full in-medium absorption (no low-px peak is in 
evidence for the critical system size in Fig. 

However, (^(Az)^^ will be significantly above 
( Az)^ ) expectation. 

Fig. 13 (right panel) also shows (z)^„^^, and 

/ d~Au 

While a ^ 30% shift is observed between (z)^_^ and 
(expected and understood through the Cronin 



effect [ig|), ((Az) 



(Az) ) for allpTmin 

d—Au \ I Au-Au 

triggers under consideration. This small system size de- 
pendence (noted previously underscores the poten- 



tial of ^(Az) j as an experimental probe for the mecha- 
nism of energy loss. 

We have verified that these results are independent 
of the inclusion of {kx) broadening within HIJING, 
but depend significantly on the inclusion of mini-jets. 
Turning off the mini-jets results in a considerable de- 
crease of ((Az)^) , resulting in a '--^ 50% difference 

\ / d-Au 

w.r.t. ((Az)^\ , as well as an expected reduction 

\ / Au-Au 

of the correlated jet sample. The divergence between 
(Az)^) and ((Az)^) in this regime is not 

/ d—Au \ I Au-Au 

surprising, since within HIJING the dominant contribut- 



ing factor to ^(Az) y is the generation of mini-jets by 
the propagating jet. We have also calculated and 



(Az)^ 



and the result is predictably indistinguish- 



In conclusion, we have illustrated the qualitative dif- 
ference between GLV and BDMPS/AMY energy loss dy- 
namics, the first being based on the decrease of the initial 
jet momentum peak, while the second predicting a broad- 
ening of the peak. We have argued that this difference 
should result in different scaling, with system size (cen- 
trality or nucleus type), of the mean and variance of the 



variable 



_ paway / pnear 



/P^ 



We have shown that in a 



able up to statistical error from d-Au. 



hadronic enviroenment, as simulated by HIJING, a sen- 
sible combination of px trigger and jet-cone angle will 
give a p°'"''J'/p^ear independent of jet-cone angle and 
momentum trigger systematics and close to the "thin 
plasma" (GLV) limit, thus providing a reliable way to 
search for the onset of BDMPS/AMY broadening. We 
have also argued that the dependence of this broaden- 
ing on system size can be related to the initial partonic 
density of the fireball. We hope that experimental mea- 
surements of (z) and ((Az)^) will be forthcoming, and 
will help us learn more about the mechanism of hadronic 
energy loss. 
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